Coherent and incoherent pumping of electrons in double quantum dots 
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We propose a new mode of operation of an electron pump consisting of two weakly coupled 
' quantum dots connected to reservoirs. An electron can be transferred within the device at zero 

\ bias voltage when it is subjected to electromagnetic radiation, thereby exciting the double dot. The 

excited state can decay by transferring charge from one lead and to the other lead in one direction. 
Depending on the energies of the intermediate states in the pumping cycle, which are controlled by 
the gate voltages, this transport is either incoherent via well-known sequential tunneling processes, 
or coherent via a inelastic co-tunneling process. The latter novel mode of operation is possible 
1/^ ' only when interdot Coulomb charging is important. The D.C. transport through the system can be 

controlled by the frequency of the applied radiation. We concentrate on the resonant case, when the 
frequency matches the energy difference for exciting an electron from one dot into the other. The 
resonant peaks in the pumping current should be experimentally observable. We have developed 
^ ■ a density matrix approach which describes the dynamics of the system on timescales much larger 

' than the period of the applied irradiation. In contrast to previous works we additionally consider 

the case of slow modulation of the irradiation amplitude. Harmonic modulation produces additional 
sidepeaks in the photoresponse, and pulsed modulation can be used to resolve the Rabi frequency 
in the time-averaged current. 



I. INTRODUCTION 



. In recent years quantum dots have attracted great attention. A quantum dot can be thought of as an artificial 
^1 ' atom with adjustable parameters. It is of more than fundamental interest to study its properties under various 
circumstances, e.g. by transport experiments By considering a double-quantum-dot system, the analogy with 
real atoms can be stretched to include artificial molecules. The analogue of the covalent bond is then formed by 
an electron which coherently tunnels back and forth between the two dots. By applying electromagnetic radiation 
' with a frequency equal to the level detuning in the double-dot system, an electron can undergo so-called spatial Rabi 
Q.^ , oscillations even when the tunneling matrix element between the dots is small. 

■ Recently, several time-dependent transport measurements on quantum-dot systems have been reported |^,^ , most 
T-H . of them being of a spectroscopic nature. It has also been suggested to construct devices from quantum dots. Examples 
\l ■ of such applications are pumps that transfer electrons one by one at zero bias voltage by using time-dependent voltages 

, to raise and lower tunnel barrier heights [0 , or systems in which coupled quantum dots (or quantum wells) are used for 
quantum-scale information processing |g[7 Several theoretical models for time-dependent transport through a double 
quantum dot have already been proposed. For instance, in Ref. Q D. C. transport was considered for arbitrary bias 

■ voltage when the signal couples to the gate voltages of the dots. At zero bias voltage the system operates as an 
. electron pump. In Ref. the D. C. current, controlled by external irradiation, was considered for finite bias voltage. 

These results were recently generalized to include time-dependent gate and bias voltages and tunnel barriers . In 
all these works a tunneling-Hamiltonian approach was used to incorporate the effects of Coulomb interaction between 
electrons on the same and on different dots. It is assumed that the barriers separating the leads and the dots are high 
^ , and therefore the wave functions have only a small overlap. A different approach would be to use scattering states of 
electrons which extend through the leads and dots which is appropriate for almost transparent barriers (see e.g. pof). 
However, the effects of Coulomb interaction are not easily taken into account in this approach. Because the transport 
mechanism in the above mentioned double-dot pumps is determined by sequential tunneling, electrons are pumped 



^ . incoherently i.e. the tunneling of an electron into and out of the device are independent events. Also, in these works 
- - ' the time-dependent signal is taken to be a monochromatic. 

In this paper we describe, firstly, an new mode of operation of such an electron pump. In this case, electrons can 
be transferred coherently through the system by means of a co-tunneling process. Our device has to be designed in 
such a way that interdot Coulomb repulsion is important. By appropriately adjusting the gate voltages the device 
can be made to operate in the co-tunneling regime or the sequential tunneling regime. The latter regime is considered 
here for comparison with previous works and should be distinguished from the coherent one. In the co-tunneling 
regime the interdot attraction between an extra electron excited by the A. C. field into one dot and the hole it left in 
the other dot stabilizes the excited state. Electrons are prevented from entering or leaving the device independently 
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and inelastic co-tunneling is the lowest non- vanishing order process for transport. This involves correlated tunneling 
events through the different junctions connecting the dots weakly to the leads. The system switches coherently from 
the excited state to the ground state each time an electron is pumped through the dots. Alternatively, this process 
can be seen as coherent transport through a double-quantum-dot qubit. The second main result of this paper is that 
we predict the effect of slow modulation of the irradiation amplitude on the pumping current. The timescale of the 
amplitude modulation is assumed to be much larger than the one set by the frequency of the unmodulated signal. 
The paper is organized as follows: in Sec. || we introduce the system. In Sec. Ill we develop the density matrix 



approach to describe the system only on a timescale much larger than the period of the applied irradiation. This 
development is similar to that in quantum-optics for resonance fluorescence e.g. and is central to our treatment of 
the problems. We apply this approach to the transport through the double-dot in the sequential and the co-tunneling 



regime. In Sec. [V we consider the case where the amplitude of the applied radiation is modulated on the large 
timescale. Finally, a summary and conclusions are presented in Sec. 

II. DOUBLE-DOT ELECTRON PUMP 

The system we consider consists of two coherently coupled quantum dots 1 and 2 connected by tunnel barriers 
to large reservoirs L and R, as depicted in Fig. |^. The leads are assumed to have a continuous electronic energy 
spectrum. The fixed difference between the electro-chemical potentials of the leads /ii = /i — eVt^ and fXR = /i — cVlj 
(e > 0) and the temperature are the smallest two energy scales of the problem. We can thus take them to be zero. 
Each quantum dot i = 1, 2 contains some number Ni electrons and is assumed to be in the ground state. We will 
concentrate on transitions between the ground states of the individual dots with different numbers of electrons. 

Disregarding all tunneling for the moment, let us consider the energy of the double dot within the standard Coulomb 
blockade model. The capacitive coupling between dot i = 1 . 2 and the attached electrodes is taken into account by 
the gate capacitance CG^and the lead capacitance c^. . The mutual capacitive coupling between the dots is described 
by the interdot capacitance ci2 . The total energy of the double dot system, when dot 1 and dot 2 are respectively in 
the Ni and N2 electron ground state, reads (see e.g. Ref. [^): 



EnuN2 = I] + ^uuN, {N, -1)+ U12N1N2 (1) 



i=l,2 1 = 1 1=1,2 

In the first term, en is the Ith effective single particle energy, ^ = 1, . . . , iV^ in dot i — 1,2: 

This incorporates the bare single particle energy e^;, a conveniently chosen offset Uii/2 and the linear shift with the 



electrode voltages. The coefficients of the voltages aic. = C^,^CGi,C(iLi = Q/cl- depend on the dot-electrode 



capacitances and the inverse capacitance matrix elements 

_ (ciC2)/q 1 _ 1 _ C12 

~ CiC2-cl,' - C,C2-C^' 

where c, = CGi + cl^ + C12 is the total capacitance of dot i = 1,2. By appropriately changing the gate voltages the 
effective single particle levels in dot 1 and 2 can be independently be shifted with respect to each other. In the second 
and third term in equation (Q), un — e'^C^^ is the intradot charging energy of dot i = 1, 2 and U12 = e^Cj^^ < un, U22 
the mterdot charging energy. 

Let us now consider the stability of a ground state of the double dot with respect to the tunneling between the 
individual dots and the leads. Assume that the gate voltages Vgi , VG2 such that in the stable state of the double 
dot there are Ni and N2 electrons in, respectively, dot 1 and 2 and denote this state by |0, 0). The stability diagram of 
the double dot near the region were this state is stable is sketched in Fig. ^ for the typical case where interdot charging 
is important: U12 < un, 1*22- The region of stability for state |0, 0) has a hexagonal shape. The stable states |ni, 712) in 
the six neighboring regions have A^i + ni and N2 + n2 electrons in dot 1 and 2, respectively, where ni ^ n2 = 0, ±1 and 
\ni + n2\ < 1. The energies of these states, En-^ n^,ni = 0, 1, are found from the right-hand-side of (|l|) by replacing 
Ni Ni + Tii. The hexagonal region is bounded by six stability constraints for state |0, 0). The first four constraints 
follow from the requirement that the energy barrier for injection {+) or emission (— ) of an electron to or from either 
lead, 
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= Ei^o - - Eq^o = eiNi+i + uiiNi + U12N2 - ^J,, 

A'^ = Eq^i - - Eofi = £2W2 + 1 + "22^2 + "12^^! - A^, 

= -1 + ^I- EoM = - {£2N2 + U22 {N2 - 1) + U12N1) , 

Al = E^i^o + M - Eq^o = (eiwi + uii {Ni - 1) + U12N2) . 



should be positive. Sufficiently far away from the four lines A J ^ = i.e. A J ^ 3> T^^/j the sequential tunneling of 
single electrons through the junctions connecting a dot and lead is suppressed, cf. Fig. ^. Here the typical tunnel 
rate is F/,/} = '^iTVL^fi\t]^ ji^ where vl^r is the density of states in the left and right lead respectively and t^ji is 
the matrix element between the states in the lead and in the dot, which depends only weakly on the energy. Two of 
sequential tunneling barriers (0) can be independently tuned by the gate voltages, the other two are related to these 



A^ = uii 



5i and A J + A^ = U22 + 82 where 5i = eiNi+i - £iN, 



-iNi + l 



e^j^, is the spacing between the 



two "active" single particle levels in dot i = 1,2. It will be convenient from here on to consider < A£ < un + Si 
and < A^ < U22 + <^2 as independent variables instead of the two gate voltages, cf. Fig. ||. Two additional stability 
constraints follow from the requirement that the energy barriers for polarizing the double dot with respect to state 
|0,0), 



£0 
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Ei.-i 



Eofl 
Eo,o 



A7 + Ai 



U12, 
U12, 



(3a) 
(3b) 



should be positive. Here A J, A^ are the positive Coulomb energies we must pay to change the number of electrons on 
each dot from the stable configuration |0,0), and — M12 is the energy we gain by creating an attracting electron-hole 
pair with respect to the stable state |0, 0). Note that the former energy also depends on uu i.e. it incorporates 
the interaction between the extra electron or hole and the electrons in both dots. Sufficiently far from the lines 
£0 = and £q = 0, cf. Fig. ||, the polarization of the double dot by a coherent co-tunnehng process is suppressed: 
£o,£o ^ ^'ct where F^^ ^ ^l.r is some typical rate for this process. In such a second order process an electron is 
injected into one dot and another electron is emitted from the other dot, effectively transporting one charge across 
the double dot. From the relation £0 + £q = X]i=i 2 ~ ^12 + Si) > we find that £o,£o > corresponds to 



Ui2 <A^+A+< E.=l,2 



U12 



ui2 in the stability diagram. 



Now consider the coherent tunneling of electrons between the dot 1 and 2. If the co-tunneling barriers are also 
larger than the matrix element T for this process, £o,£o 3> T, then the polarization of the double dot by coherent 
tunneling of an electron between the dots is also suppressed. Under these conditions the D.C. transport through the 
double dot is blocked at low bias voltage i.e. we have the Coulomb blockade. 

This situation is changed, however, if we apply electromagnetic radiation to the system. Assume that a time- 
dependent oscillating signal is present on the gate electrodes which will shift the energies of single levels without 
altering the wave functions too much, so that the time-dependent energy difference between states |0, 0) and | — 1, 1) 
becomes 



e{t) = £0 + V cos ujt, 



(4) 



where V is the amplitude and lu the frequency of the externally applied signal. When the frequency of this applied 
radiation matches the time-independent energy difference £0 between states |0, 0) and | — 1, 1), it is possible for an 
electron from the left dot to tunnel to the right one by absorbing one energy quantum lu Sq from the field. In 
principle, this electron can now leave the system by tunneling to the right lead, resulting in the state | — 1,0). An 
electron from the left lead can then tunnel to the left dot, thus restoring the ground state. Effectively, an electron has 
now been transferred from the left electrode to the right one. Alternatively, the electron can coherently tunnel back 
by emitting an energy quantum resulting in state |0, 0). This transport cycle, |0, 0) ^ | — 1, 1) | — 1,0) |0,0), 
is not the only one. Another possible sequence, in which the system passes the intermediate state |0, 1), is given by 
|0, 0) I — 1, 1) ^ |0, 1) — > |0, 0). The tree other states can be disregarded under the following conditions. Firstly, 



the field should not be resonant with the transition to the other excited state |1, —1): |£o 



< £n 



This is the 



case when the distance between the two excited levels is much larger than the frequency detuning of the frequency, 
|£o — £o| ^ \Soj\ = \eo — uj\, which, in the stability diagram, corresponds to (cf. equations ^ 
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1=1,2 



> \Suj\ 



(5) 
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Secondly, the states |1,0),|0,— 1) should not be resonant with intermediate states of the transport cycle |0, 1) and 
I — 1,0), respectively: li^i^o — i?o,i| > T and |i?o,-i — £^-i,o| > T corresponding to 

|uii + ,5i-(A-+A+)| »r, (6a) 
|u22 + 52-(A^+A+)| >r. (6b) 

Thus under these conditions an electron can be excited from the left dot into the right one, but the probability of 
exciting an electron from the right dot to the left dot can be disregarded. 

The details of the transport mechanism of a pumping cycle depend on the energies of the intermediate states | — 1, 0) 
and |0, 1) relative to the pumped state | — 1, 1) as shown in Fig. || which are controlled by the gate voltages. The 
energy barrier for injecting an electron from the left lead into dot 1 and for emitting an electron from the dot 2 dot 
to the right lead, 

A+ =i;o.i =ui2- A;^:, (7a) 

A^ = i;_i,o + Ai-^^-i4 = ^*i2- a;;j, (7b) 

can be positive or negative, depending on the position in the stable region of |0,0). Here uyi is the energy we must 
pay to break up the attracting electron-hole pair with respect to state |0, 0) and — A^, — A^ is the Coulomb energy we 
gain by changing the number of electrons on one of the dots to the value of the stable configuration |0, 0). In this paper 
we consider two regimes of operation of the double dot electron pump, which are schematically depicted in Fig. ||: (I) 
both barriers are negative, Aj, <^ — F^, — F^: the pumped level can decay through sequential tunneling processes; 
(II) both barriers are positive, AJ,A^ ^ Fi,,Fj^: the excited state is stable with respect to sequential tunneling 
ay through an inelastic co-tunneling process. We do not consider the more complicated intermediate case 
^ rLjT-R where resonant processes between leads and dots are important. In the stability diagram in 
wo regimes correspond to 

(I) |A--A+| < (A- + A+)-2ui2, 

2wi2 < A^ + A+ < {uu - ui2 + S,) , (8a) 

4=1,2 

(11) |A^-A+| < 2^.12 -(AZ+A+) , 

ui2 < A^ A+ < 2^12, (8b) 

where < stands for "separated by energy large compared with Tl h". The narrow strips defined by conditions (|^) 
and (^) should be excluded from these regions. 

In regime I (Aj, Aj^ < 0) the charging of the individual dots dominates the transport. The system relaxes to the 
ground state via the two sequential (and thus incoherent) tunneling processes described above. As shown in Fig. ||(a), 
the four tunnehng processes are described by the rates F)^, F)j, F^, and F^j, respectively. The following rate equations 
describe the density matrix in the sequential tunneling regime: 

dtP-io-m = +rV-ii -11 ~ ^lP-10-10 (9a) 

<9tPoo,oo = -ir(/9_ii,oo - Poo,-ii) + r° p-io -10 + ^nPoim, (9b) 

9tP-ii,-ii = +iT(/9_ii,oo - Poo.-ii) - (r?;, + ri) P-11,-11, (9c) 

dtPoim = +F^p_ii._ii - F]j/9oi,oi, (9d) 

<9tp-ii,oo = -iT(/9oo,oo - P-ii -ii) - ie (t) P-ii.oo - ^ {^i + ^r) P-n^oo, (9e) 

Here, and throughout this paper, units are used such that h — 1. The diagonal elements give the probabilities for an 
electron to be in the corresponding states and probability is conserved i.e. at any time t 

P-10,-10 + Poo, 00 + P-11,-11 + Poi,oi — 1- (10) 

The non-diagonal elements p_ii,oo — Poo,-ii describe the coherence between states | — 1,1) and |0,0). In general 
the tunnel rates F^'^ and F^^ depend on the energy difference between the states of the transition. We can take 
F^'^^ = Flu = 2i:vL ji \tL,R^ when the density of states vl.r in the left and right lead respectively and the matrix 
element t^ ji between the states in the lead and in the dots depends only weakly on the energy. The average current 
through the system is given by 
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^/e = ^rP-11-11 + ^rPoi.oi 



(11) 



In regime II (Aj,A^ > 0) the interdot attraction of the electron- hole pair is dominant over the individual 
(de)charging of the individual dots. The decay of the excited state | — 1,1) via sequential tunneling is blocked 
as shown in Fig. ||(b). However, transport is still possible via inelastic co-tunneling of electrons |l^]. When the 
system is in state | — 1, 1), two electrons can tunnel simultaneously through different barriers, one going from the left 
lead to the left dot, and one from the right dot to the right lead. Because in this transport process a state is virtually 
occupied these two tunneling events cannot be treated independently. The necessary energy is provided by relaxing 
the system to the ground state |0,0), thereby releasing an energy E « eo-The co-tunneling rate can be calculated 
with Fermi's Golden Rule. The relevant matrix element is a sum of matrix elements for the two possible coherent 
processes which transfer one electron from L to R. The co-tunnel rate is obtained by integrating the partial rates for 
transitions over the different final states of the leads which are assumed to be uncorrelated: 



Fct = 2myLVfi / dsL 

-oo J —oo 



deRf{eL){l-f{eR)) 



thtR. 



ihtR 



A+ - EL Afl, + Er 



5 (eo + £l- £r) 



(12) 
(13) 



Here the matrix elements ^l.r for tunneling through the left and right barrier and the densities of states vl,r in the 
left and right electrode are assumed to be energy independent. The zero-temperature co-tunnel rate in our electron 
pump is 



Tct = 



rr.r 



R 

27r 



£0 



A 



£0 A^ 



a; 



- £o 



ln(l + |^)+ln(l + || 




(14a) 



(14b) 



Note that from ep = Ui2 — ( Aj + ^_r) (cf- equations ( |3a| ) and (0)) we observe that within regime II we can have 



eo ^ AJ, Ajj. The energy denominators in (14a) reflect the fact that the tunneling occurs via the virtual occupation 
of two states. In contrast to the incoherent sequential tunneling mechanism, the only relevant density matrix elements 
are those between states |0, 0) and | — 1, 1) since states | — 1,0) and |0,1) are occupied only virtually. Taking the 
co-tunneling processes into account we obtain the equations of motion for the density matrix elements: 



where p. 



11,00 



9(Poo,oo = ~i2^(p-ii,oo — Poo,-ii) + TctP-ii.-ii 
i9tP-ii,-ii — +ir(p_ii,oo ~ Poo,-ii) ~ TctP-ii.-ii 
dtp^ii,m = -ir(poo,oo - P-ii.-ii) - i£ it) p-ii,oo 

Poo -11 ^^"i probability is conserved: 

Poo, 00 + P-ll.-ll ~ 1 



^rctPoo, 



■11; 



(15a) 
(15b) 
(15c) 



The current is 



/(i)/e = retp-ii,-ii(t). (16) 

In both regimes the irradiation relaxes the constraint of energy conservation during tunneling by allowing an 
electron to absorb or emit a multiple of the energy quantum uj. This gives rise to an enhancement of the zero bias 
D.C. component of the current and additionally it introduces fast oscillations with small amplitude which do not 
interest us. In the next section we show how to extract the slowly varying component of the density matrix from 
the exact equations (H) and ( p^ respectively. We point out that above we have written expressions for the particle 
current only. The displacement current can be disregarded here since it does not contribute to the D. C. current. 
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III. TIMESCALE SEPARATION 



In this section we consider the dynamics of the two coherently coupled levels \a) — |0, 0), |6) = | — 1, 1) on timescales 
much larger than the period of the applied irradiation = 2tt /uj. The details of the other states in each regime 
are only important for the incoherent processes which depend only weakly on the time-dependent energy difference 
between basis states \a) and |6) 

e{t) ^ eo + V [t) cos {ijjt) (17) 

Here we also allow the amplitude V{t) to be modulated on a timescale which is large relative to t^. The timescale 
separation for both regimes can be done in the same way. The coherent part of the dynamics of the state of the 
system depends only on the time-dependent: H (t) = Hq (t) + Ht where 

no{t)^\e{t){\a){a\^\h){h\) (18) 

introduces the energy difference between double dot states with zero extra electrons and Ht describes the tunnel 
coupling between the dots: 

nT^T{\a){b\ + \b){a\) (19) 

We assume that the tunneling amplitude is much smaller than the time- independent energy difference T Eqj whereas 
V (t) can be of arbitrary magnitude. The tunneling to and from the reservoirs brings the system into a mixed state 
which can only be described by a density operator p which obeys the Neumann-Liouville equation with dissipative 
terms added to the right-hand-side ||l^ , |l5|| : 

dtp = -i [H, p] + CincP- (20) 



Since the incoherent part in equation (^ and (y^ is invariant under a phase transformation of the nondiagonal 
elements, we can derive from (|2^) a set of equations that describes the dynamics on large timescales by first performing 
a standard time- dependent basis transformation |Q on the density matrix. A rapidly varying time-dependent phase 
factor is absorbed in the nondiagonal elements of p 

Pab = (21) 

and the diagonal elements are left unchanged. In this new basis the generalized Liouville equation for p takes the same 
form as equation ( ^o|) with the same incoherent part and a new Hamiltonian Ti' [t) ~ TYq + {t) with a renormalized 
energy difference and a time-dependent tunnel amplitude 

K = \{e{t)-dtm){W){a'\-\b'){b'\). (22) 

H't {t) = re-'^(*) (|a'> {b'\ + |6'> (a'|) . (23) 

We choose the phase to be {t) — nuit + {V{t)/uj) sin(ijj<) to obtain time-independent diagonal elements 

e{t) - dt<j)it) = £0 - nuj (24) 

which vanish at the n-photon resonance nuj = eg in which we are interested. Furthermore, on small timescales the 
new tunnel matrix element is a periodic function of time and can be expanded in a Fourier series: 

H't (t) - e-'("+")"*H^^^ (t) , (25) 

m— — oo 

Wt,„ it) = Jrn (^) T da') {b'\ + \b') (a'l) . (26) 

Likewise we expand the density operator into rapidly oscillating contributions with amplitudes which vary on large 
timescales: 

oo 

p'(t)=^p'W(t)e--*. (27) 
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Inserting this into the generalized Liouville equation for p' we obtain an infinite set of coupled equations for the 
slowly varying coefficients /5'('"\ The amplitude of the fast oscillations p'^^ is of order T jio ~ T /eq <C 1 and can be 
disregarded: the "D.C." component (t) satisfies 

= -i \h'^ + , + (28) 

Thus the nearly isolated states |a), \h) irradiated at a resonant frequency lu so/n are equivalent to almost degenerate 
states \a'), \b') coupled by a tunneling matrix element 

T it) - J-n (^) T = (-If J„ (^) T (29) 

which only varies on large timescales. In the following we will only consider the 1-photon resonance i.e. to ~ Eq 
and we omit the superscripts used above to distinguish the slowly varying components from p itself. However, the 
equations can be generalized to the n-photon case by replacing Ji (V/lj) J„ (V/oj) and Eq — uj ^ Eq — noj. Due to 
the oscillatory behavior of the Bessel function Ji the coherent tunneling amplitude can be tuned between < T < 
0.58T by varying the amplitude/frequency ratio of the irradiation over a range < ^ ^ 1.84ti;. The vanishing of the 
effective tunnel matrix element, T = for non-zero V/uj is one of the features which distinguishes photon-assisted 
tunneling from adiabatic electron transfer. A similar renormalization of the tunnel coupling to zero is also known 
from driven double well potentials |p^ . 

The approach developed above allows us the extract the slowly varying components of the current in both regimes 
of operation of the electron pump. We point out that we have only changed the coherent part of the dynamics which 
is the same in both regimes. 



A. Sequential tunneling regime 

The equations (|^) describe the dynamics of the double dot in the sequential tunneling regime on the long timescale 
when we replace e (i) ^ eo — ijJ, T ^ T = — Ji (^^-^^ T. The solution of these equations tends to a stationary value, 

which is independent of the initial conditions, on typical timescale max{(r^'^)~^, (F^^)^^}. The solution of dtp (t) — 
gives a Lorentzian lineshape of the current peak as a function of Eq ~ 

I/e^ I^,^wy[w'' + {eo-u;)% (30) 

where the maximum current /max and half-width at half-maximum w are given by 

/max/ s — 1/ I r'l+r]^ ^ ^ (rf' r^) 4^2 " ^ i (31) 



w 



^ J? I ^ t 



T2 + ^lUdpil. (32) 



Let's assume that the tunnel rates for each barrier are the same: T^'^^ = T^l.r- If the double dot is weakly coupled 
to the leads i.e. Tl.r ^ T <C £0 ~ ^ then one can access the regime Tl^h, \eo — uj\ <S^T hy adjusting the irradiation 
amplitude to ^ « 1.84^, where the coherent state in the double dot dominates the transport properties . The two 
delocalized states in the double dot are independent channels for transport and the current increases with T^ ji: 

Ira../e = 1/ ( + j^) , W^2f (33) 

In the opposite regime |eo ~ "^l T ^l,r, which can be accessed by tuning V ^ 1.84a), the decoherence due to 
tunneling to and from the reservoir dominates the transport. In this case the height of the current peak is proportional 
to and decreases with enhanced tunneling Tj^ ji: 

/max/e = 4rV {Tl +Tr), w = (T l + Tr) /2 (34) 

The current peak /max/e reaches a maximum T/2 as a function of T^ jf when = Tr = 2T. This can be understood 
as the precise matching of tunneling times: the time for half a Rabi oscillation in the double dot is exactly equal to 
the time for filling the left dot and for emptying the right dot. In Fig. ^(a) we have plotted the maximum current and 
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the width as a function of the tunnel rate F relative to the tunnel coupling T. If the double dot is strongly coupled 
to the leads i.e. T <^ Tl,r <C Eq ~ only the regime T, jeg — w| <C Tl,r is accessible. 

We point out that in the sequential tunneling the transport on "large" timescales is equivalent to transport of "free" 
electrons (i.e. negligible interdot repulsion) through a double dot with renormalized static parameters Eq Eq — uj, 
T —f Ji (V/uj) T, h Also, the results here are similar to the analytical results obtained for the non-interacting 
case in [||, where only sequential tunneling. 



B. Co-tunneling regime 

In the co-tunneling regime two well-separated timescales are involved, namely the "long" time tct = between 
two co-tunneling processes and the "short" time of the process itself ivirt = £(7^ during which the energy is uncertain 
by an amount Eq. Near resonance the applied frequency matches the detuning of the levels so tvirt ~ 27r/w = t^. 
The density matrix approach developed above describes the system on the "large" timescale t^- The equations (|l5| ) 
describe the dynamics of the double dot in the sequential tunneling regime on the long timescale when we replace 
e {t) ^ Eq — U!,T ^ T = — Ji ^^-^^ T. The stationary solution of these equations gives a Lorentzian current in Eq — lu 
with height /max and half-width w: 

^max/e = 1/ (t^ + i#) , (35) 

w = ^2f2-H%. (36) 

For weak coupling to the reservoirs Fct ^ T the peak height is constant and the width increases linearly with the 
coherent coupling T: 

/max/e = ^Fct, W = V2f 

In the opposite case Fct ^ T the peak height is small but increases rapidly with T whereas now the width is constant: 

/max/e = 4f'^/Fct, W = ^Fct 

The current peak /max/e reaches a maximum T/V2 as a function of Fct when Fct = 2V2T. In Fig. |(b) the scaled 
pumping current is plotted for different values of Fct /T. In the co-tunneling regime electrons are transferred through 
strongly correlated transport channels. Therefore the condition for the maximum current cannot be understood as a 
the precise matching of tunneling times as in the sequential tunneling regime (factor \/2). Equations (^5|) coincide 
with those that describe the transport of electrons which are correlated by strong Coulomb repulsion through a double 
dot at high voltage bias, with renormalized static parameters eo — > eo — w, T ^ Ji {V/oj) T, Fct — > F/j in the limit 
where the tunneling to the into the left dot is so fast, F^ ^ F/j, T, that the current doesn't depend on it anymore 0. 
In this system the correlation of the conduction channels prevents more than one channel from being occupied and 
reduces the effective tunnel rate by a factor 2 in ( p5| ) Jl7[ . This blockade can be obtained formally from the density 
matrix equations of the "free" electrons (^ by taking F]^ ^ — > and F° 3> F^j = Fct • 

Comparing the co-tunneling and sequential tunneling regime the main difference is that for fixed F^ /j the co-tunnel 
rate oc F^F/j is much smaller: Fct ^ ^l,r- For F^^^j ~ T the sequential tunneling current can be near its maximal 
value ~ 2T , whereas the co-tunneling current will be cx Fct ^ T. However, for Tj^^n ^ T the sequential tunneling 
current is much smaller than T and it is possible to adjust Fct ^ T so that the co-tunneling current takes its maximal 
value ~ 2\/2T which is much larger. Thus for fixed T the maximal co-tunneling current is larger by a factor 
than the maximal sequential tunneling current (Fig. ^ . The width of the co-tunnel peak is also smaller than in the 
sequential tunneling regime. 



IV. MODULATED IRRADIATION OF A DOUBLE DOT 



In this section we consider the sequential tunneling regime already discussed in section III A and apply our approach 
for large timescales to the cases of pulsed and slow sinusoidal modulation of the irradiation amplitude. 
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A. Response to irradiation pulses 



By means of irradiation pulses quantum states in the dots can be manipulated. Assuming the double dot to be in 
the ground state 



Poo, 00 — liP-10,-10 — Poi,oi — — P-iim — Poo, -11 — 



(37) 



at t = we solve equations with e (t) Eq — T — > T = — Ji ( ) T for the time-evolution under the influence 



of the irradiation for the case r°'^^ = F: 



I{t) = I 



1 - e 



-r* 



cos i^Rt) + — — sin {flnt) 



(38) 



where f^R = y {eq — uj) + AT"^ is the Rabi frequency. For T ^ the solution tends to the stationary current (|3 
derived before: 



//e = l//'iilL + ^ 



(39) 



When the irradiation is switched off at t = Tp, the current decays exponentially as e"'"'*^'^*'-' to zero, as can also be 
seen by solving the equations for the case lo = T ~ {). One can resolve the Rabi oscillations in the D.C. -current by 
considering the current averaged over a series of identical pulses p8| with delay r^, as a function of the pulse- length 
Tp (see inset of Fig. ||): 



^dc (Tp) = — , , 



I it) dt + I (Tp) 



(40) 



Here — Tp ^ F ^ to ensure that the system is prepared in the ground state (|3^) at the beginning of the pulse. In 
the physically interesting case of weak coupling to the leads F ^ £7^ we obtain 



Idc (Tp) « I — 
Td 



1 — e ^■^f cos (r^RTp) 



(41) 



The oscillations are most clearly resolvable for Tp ^ F ^ (Fig. ||). The period of the coherent oscillation at resonance 
u! = £oi 27r/fiR = tt/ ( Ji {Vq/ui) T), can be tuned by varying the irradiation power Vq. 



B. Sinusoidal amplitude modulation 

Now consider the case where the amplitude of the irradiation is slowly sinusoidally modulated with small modulation 
amplitude V <^Vo: 

V{t) ^Vo + Vcos{nt). (42) 

The case where the modulation amplitude is of the order of or larger than the irradiation amplitude is physically not 
very interesting because the system then exhibits trivial Fourier peaks at lu, lo + Q and uj — Q as a function of Eq in 
the D.C. -current. To find an analytical solution, wc rewrite equations (||) with e (t) ^ Sq — uj,T ^ T — —Ji (^^-^^ T 
in matrix notation: 

%-{t^f)p^c. (43) 

where p = (p-io.-io,Poo,oo, P-ii,-ii, Poi,oi, P-ii,oo)^, <? = (F)j, 0, 0, 0, 0)^and expand the Bessel function Ji {V (<) /w) 
in a Taylor series to second order in V . If we now consider the Fourier coefficients p„ of p{t) and r„ of Ji {V {t) /oj) T, 
we find the following equations for the D.C. -component and first harmonic, if we disregard higher harmonics 
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0= (f + Tof) po+Tif (p+1 +/7_i) + c, (44) 
p±i = -TJt + TofTi^l) 'fpo, (45) 



from which we can solve for the D.C. -component po- If we furthermore assume that T£ = F <C T, the solution 
shows additional sidepeaks in the photoresponse of the system, which are in good approximation Lorentzians along 
the hyperbola Eq — uj = Vfl'^ — AT"^ i.e. — SIr as plotted in Fig. (^. The height of these peaks is 

I^T[^^l^^^^^^^ (46) 
with a = n/T,T ^ Ji{Vo/u!)T and the half- width at half maximum is 




(47) 

These sidepeaks thus resolve the Rabi splitting as described in Refs. ||6|,|l9| in terms of quasi-energies. The height of 
these peaks can be of the order of the first satellite peak and should therefore be experimentally observable. Notice 
that ai Eq — = 2T the width diverges; however, the current peak at this point vanishes since the exact matching 
Eq — uj gives a resonance to which the modulation cannot add extra current. 

Thus by applying a high frequency lu we reduce the energy spacing by a large amount, Eq Eq — to, and modify 
the tunneling matrix element with an intensity {V) dependent factor, T ^ T. The lower frequency allows one 
to precisely match the remaining small energy difference Eq — uj without significantly altering the tunneling matrix 
element T, thereby inducing a photo-current. 



V. CONCLUSIONS 



We have considered an electron pump consisting of an double quantum dot subject to irradiation. An incoherent 
and a new coherent pumping mechanism were discussed. We have derived equations of motion for the density matrix 
elements of the double-dot system which are time-averaged over an interval which is long compared to the period 
of the applied signal. From these equations we calculated the pumping current in both regimes. In both cases the 
current peak is a Lorentzian. Surprisingly, for fixed effective tunnel coupling the maximal pumping current in the 
co-tunneling regime is larger by a factor \/2 compared to the value in the sequential tunneling regime, where the 
maximum occurs at a different value of the tunnel rates for each regime. Experimental realization of this device 
would allow for a systematic study of coherent transport through a solid-state qubit. Moreover, modulation of the 
irradiation amplitude exhibits interesting phenomena: a train of short pulses should resolve the Rabi frequency in the 
time-averaged current as a function of the pulse length, and sinusoidal amplitude modulation should provide a tool 
to resolve the Rabi splitting. 
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FIG. 1. Double dot electron pump. The asymmetric gate voltages allow one to induce a photo-current at zero bias. One can 
excite an electron from the highest level filled in the left dot to the first empty level of the right dot with an AC field that is 
resonant with this transition, but not with the corresponding transition in the opposite direction. 



FIG. 2. Stability diagram of the double dot system near the region where the state [0, 0) is stable, with A'^i and A'2 electrons 
in dot 1 and 2, respectively. Along the horizontal and vertical axis the sequential tunneling barrier and A J, respectively, 
are varied independently. 



FIG. 3. Energy diagrams of the double dot electron pump together with the leads: a) sequential tunneling regime, b) 
cotunneling regime. The Coulomb interaction between the extra positive charge on the left dot and negative charge on the 
right dot determines the position of level £-1,1 relative to -B-1,0 and -Eo,i- 



FIG. 4. Lorentzian current peak height (bold solid line) and width (dashed line) as a function of a) F/T (F = Fl = Tr) in 
the sequential tunneling regime and of b) Fct/T in the cotunneling regime. For fixed T the maximal co-tunneling current as a 
function of Fct is larger than the maximal sequential tunneling current as a function of F. 



FIG. 5. Double dot subjected to an irradiation pulse train. Figure: Current averaged over a series of identical pulses,in units 
of number of electrons transferred per pulse, as a function of the pulse length for eo — lu = 30, T = 1, F = 0.1 and several 
irradiation amplitudes V. Inset: time dependent current during one pulse. 



FIG. 6. Electron pump with harmonically amplitude-modulated irradiation. Current as a function of the renormalized level 
detuning (resonance mismatch of basis frequency) and the frequency of the amplitude modulation for T = 1, Fl.r = 0.2, 
Vo = 30, V = 10. 
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This figure "sinusmod.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0004139v2 
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